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Abstract. We show that under certain simple assumptions on the topology (structure) of networks 
of strongly interacting chaotic elements a phenomenon of long range action takes place, namely that the 
asymptotic (as time goes to infinity) dynamics of an arbitrary large network is completely determined 
by its boundary conditions. This phenomenon takes place under very mild and robust assumptions 
on local dynamics with short range interactions. However, we show that it is unstable with respect to 
arbitrarily weak local random perturbations. 

1 Introduction 

Despite a number of impressive results in the field of coupled map lattices (CML) up to now only the 
situation when the interaction between individual systems is very weak is reasonably well understood 
from the mathematical point of view (see a review of recent results, e.g. in [3J Oj]). In that case 
it is known that under mild technical assumptions the asymptotic behavior of a CML is very close 
to the behavior of a direct product of individual uncoupled systems. Note that counterexamples 
demonstrating very different types of behavior under arbitrary weak interactions are known as well 
(see J3j for details). On the other hand, numerical studies (see, e.g. (JJ) suggest that when the 
interactions becomes stronger a kind of phase transition from one to several invariant SBR measures 
takes place. Despite a number of efforts during last ten years to give a mathematical description of this 
phenomenon no rigorous results for infinite systems have been obtained. Probably one of the reasons 
is that all mathematical results in the field are of perturbative nature and corresponding methods do 
not work well for large perturbations. Another problem here is the necessity to work with infinite 
collections of individual maps. Therefore, it is harder to study corresponding functional spaces. 

In this paper we shall follow a more "physical" approach to the analysis of CMLs' dynamics in a 
strong interactions case. Indeed, a general receipt in statistical physics to an analysis of an infinitely 
large system is to consider its restriction to a finite volume (box) with certain boundary conditions and 
to study what happens when the size of the box grows to infinity. The case when the dynamics inside 
of the (sufficiently large) box is not influenced by boundary conditions corresponds to the situation 
known as the case (range) of weak interactions (or as the absence of phase transitions). We shall be 
mainly interested in the opposite situation when the asymptotic (as time goes to infinity) dynamics 
inside of the box (with an arbitrary large diameter) is completely determined by the boundary con- 
ditions and which is reasonably to call the long range action (LRA). As we shall see the LRA takes 
place under sufficiently strong interactions in networks whose topologies satisfy the property which 
we call unidirectional and which generalizes well known one-dimensional CMLs with unidirectional 
interactions |Hj. 

The paper is organized as follows. In the next section we give some technical definitions and formu- 
late our main result about the existence of the LRA in a reasonably large class of CMLs (Theorem l2.1|) . 
Sections |3] and ^ contain the proof of this result and the proof of the equivalence of the notion of the 
unidirectional interaction to the absence of cycles in the directed connectivity graph of interactions in 

'Russian Academy of Sciences, Inst, for Information Transmission Problems, and Observatoire de la Cote d'Azur, 
e-mail: blank@obs-nice.fr 

' School of Mathematics, Georgia Institute of Technology; e-mail: bunimovh@math.gatech.edu 



1 



a network. Section|5]is dedicated to the discussion of the topology of networks corresponding to CMLs 
under study and to the influence of different types of boundary conditions on the existence of LRA. 
In Section El we study random perturbations of local dynamics and show that the LRA phenomenon 
is stochastically unstable. Finally in Section [7] we briefly address some open questions and compare 
our results to some situations known in the theory of synchronization. 



2 The setup 

Let T : X — > X be a map from a compact metric space (X,p) into itself satisfying the Lipschitz 
condition with the contracting constant At, i-e. 

A T := sup P -S^Iyl < oo, (2.1) 
x^y&X P\x, y) 

(e.g. an Anosov map on a finite-dimensional torus with the Euclidean metric). Then the pair (T, X) 
defines a dynamical system. Consider an infinite collection of dynamical systems (Tj, X) with i € Z rf 
and a certain interaction between them. We always assume that 

Ay. < At < oo 

for all i G Z d and denote the direct product of the systems (Tj, X) by (T, X). The maps T\ we shall call 
local maps. Elements x of the infinite dimensional space X will be called vectors and their coordinates 
will be denoted as Xi := (x)i £ X with i £ 

Definition 2.1 By the interaction on the lattice X := X^ d we shall mean the map X : X — > X 
preserving the 'diagonal', i.e. 2(Ce) = Ce for any constant C and the vector e£l with all entries 
equal to 1. We restrict ourselves to the case of continuous interactions satisfying the assumption that 
there exists a constant A^ < 1 such that 



p((lx) e , (iy)e) < A f Yl PfaVi) ( 2 - 2 ) 



for any £ £ Z d and any vectors x, y € X. 

In particular, this means that for any pair of vectors which differ at only one coordinate the map 
T is a pure contraction. It is tempting to assume that the interaction I is a pure contraction in the 
product metric, which would simplify a lot the further analysis. At the first sight it seems that the 
difference between the assumption (|2.2|) and the pure contraction is not essential. However a closer 
look shows that the latter contradicts to the definition of the interaction which says that the map 2 
preserves the 'diagonal', i.e. the distance between the vectors X{C\e) and X{C2&) is the same as the 
distance between the vectors C\e and C2& for any constants Ci,C2- 

The strength of the interactions we shall measure by the difference 1 — A^. Indeed, without the 
interactions one has = 1. A "typical" example: I is a linear map represented by an infinite 
matrix with normalized rows. Local interactions correspond to the case when only terms near the 
diagonal of the matrix differ from zero. Again it is useful to have in mind a "standard" example of 
the local interaction (often called a "diffusive coupling") on a one-dimensional lattice when after the 
interactions the new value at the site i is calculated according to the formula: 

Xi -> (1 - 2c)xi + c(xj_i + x i+ i). (2.3) 

Here the interaction is parametrized by the value < c < 1/2 and A^ = max(l — 2c, c). 

Another important and well known example describes unidirectional interactions on a one-dimensional 
lattice defined by the formula: 

Xi -> (1 - c)xi -\-cxi-i. (2.4) 
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In this case < c < 1 and = max(l — c, c). 

Both these examples belong to the class of short range interactions when subsystems interact 
only if the Eucledian distance between the corresponding sites do not exceed a certain threshold. To 
emphasize that our results do not come as a consequence of long range interactions (when subsystems 
being arbitrary far apart may interact with each other) we shall restrict the consideration only to 
short range interactions. 

Definition 2.2 The superposition T := To T of the interaction and the direct product of local maps 
defines the dynamical system (T, X) which is called a coupled map lattice (CML). 

Definition 2.3 A box B is a bounded subset of the integer lattice 7L d . 

To restrict the dynamics to the box B we define it as follows: 

^*).:={(( f0T >), lf!£S (2.5) 
' 1 yxi otherwise. 

Definition 2.4 A collection of values of x outside of the box B is called boundary conditions and is 
denoted by where by B c we denote the complement to the set B. 

If we choose certain boundary conditions the map 7jg defines a finite dimensional dynamics inside 

a finite box B. The collection of functions Xi{t) := ^7jg ^(0)^ with i € B and t E Z + will be called a 

solution of our problem corresponding to the initial data x(0). Clearly x(0) can be decomposed into 
the genuine initial data x^(0) and the boundary conditions x^c. 

Definition 2.5 By the long range action (LRA) we shall mean that for a given box B and given 
boundary conditions X|g c there exists a vector x\jg (called a limit solution) such that for any initial 
data x\g(0) the corresponding solution x(t) converges pointwise with time t to x\q. 

The main issue here is that the boundary conditions define the long term dynamics inside the box 
B in a unique way contrary to the usual situation when the dynamics inside the box is asymptotically 
independent on the boundary conditions. On the other hand, we do not assume that different boundary 
conditions lead to different solutions. The point is that if the local maps are not bijective a pair of 
very different boundary conditions may lead to the same limit solution. On the other hand, under 
some simple additional assumptions (of non-degeneracy type) close but different boundary conditions 
lead to different limit solutions. 1 

Definition 2.6 We shall say that the interaction Z in the box B is unidirectional if there is a positive 
function M : B — > Z 1 such that 

- Mi ^ j € B we have M{i) ^ M(j); 

- Vx € X and Vi € B the value (lx)i does not depend on the values Xj for all j with Af(j) > M(i). 

In other words the function N is an enumeration according to which sites in B with larger numbers 
do not 'interact' with sites with smaller numbers. 

It is convenient to describe a topology of interactions in terms of connectivity graphs. 

Definition 2.7 A connectivity graph G := G(B,I) of the interaction Z in the box B is a directed 
graph G, whose vertices correspond to the sites of B. Besides the edge (i,j) connecting the vertex 
i € G to the vertex j € G belongs to G if and only if the value (Zx)j depends on the value Xi. 
The latter means that the value (Zx)j assumes at least two different values when we are changing Xi 
assuming that all other coordinates of x are fixed. 

x See a discussion in the end of Section [2] 
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Figure 1: Examples of unidirectional connectivity graphs 

To take into account the boundary conditions we include to the graph G all edges leading from B c 
to the points belonging to the box B and extend the function TV to these sites setting its value to 
on them. 

For each site £ £ B we denote by T(£) the collection of sites in the box B from which there are 
directed paths in the graph G leading to £, and denote by L(£) the length of the shortest path in G 
leading to £ from the points outside the box B. If there are no such paths we set L{£) = co. 

Now we are ready to formulate the main result of the paper. 

Theorem 2.1 Let the interaction in a finite box B be unidirectional and let for any i G B there be 
a path in G(B,Z) from a site outside of the box B to i. Then for any boundary conditions the LRA 
takes place whenever A^Ay < 1 . 

Note the local maps may be both regular and chaotic here. The simplest example of the unidirec- 
tional interaction is a linear one-dimensional (i.e. d = 1) system of unidirectionally interacting maps 
with the interaction described by the relation ()2.4j) . Of course, if this situation would exhaust the 
list of examples there would be no reason to discuss them at all. It turned out that there is a much 
bigger class of systems with effectively unidirectional interactions. This class of interactions have a 
very simple and intuitive interpretation in terms of a connectivity graph. 

Theorem 2.2 An interaction is unidirectional if and only if the corresponding directed connectivity 
graph has no cycles. 

Two examples of graphs satisfying this property are depicted in Fig.^ The left panel in Fig. ^cor- 
responds to the interaction similar to the one which is well known under the name Toom's North-East 
voting model in statistical physics (see e.g. 0), while the right one demonstrates a bit more complex 
topology of interactions. Observe that in both examples the graphs may be continued periodically to 
an arbitrary large rectangular box. 

3 Proof of Theorem EH] 

First we prove the following technical result. 

Lemma 3.1 Let {vi(t)}^L 1 be a collection of functions Vi(t) : Z^ —* X with the parameter t € Z^ and 
let the function u(t) : l} + — > X be defined by the following (non- autonomous) relation: 

U(t+1) :=$(T u(t),T 1 V 1 (t),T2V 2 (t),...), (3.1) 



4 



where u(0) := no £ X, the maps r\ : X — > X satisfy the inequality A n < At for all i, and the map 
<3? : X N+1 — > X satisfy the inequality 

N 

p($(z , *(4 4, •••))< A* (3.2) 

i=0 

/or any too sequences {zi}, {z^}. Then whenever A$At < 1 f/iere exists a function u(t) : — > X 
(depending on V{) such that for any uq G X the solution u{t) of the equation \3. 1\) with u(0) = Uo 
converges with time at exponential rate to u(t), i.e. p(u(t),u(t)) < C^ 1 for all t > and some 
7 < 1, C < oo. Moreover, if there exists a vector v € X N such that maxi<j<7v p(vi(t),Vi) < Crf for 
all t > and some 7 < 1, C < 00, then the limit u := lim^ 00 u(t) exists as well. 

Proof. Choose two different initial values uq , u' Q G X and denote the corresponding solutions of the 
equation (|3.1jl by u(t) and u'(t) respectively. Then 

p(u(t + l),u'(t + l)) 

= p($(T u(t),T 1 V 1 (t),T2V 2 (t),...), ^(T u'(t),T 1 V 1 (t),T 2 V 2 (t),...)) 

< A^p{T u(t),T u'(t)) < A^A T p(u(t),u'(t)) < ... < (A»A T )V(«o>«d) ^ °> 

which proves the first assertion with C := max x ^ e x p(x,y) and 7 := A$Aj- < 1. 
To prove the second assertion consider 

p(u(t+l),u(t)) = p($(T u(t),T l V l (t),T2V 2 (t),...), ®(T u(t- l),T lVl (t- l),T 2 V 2 (t- 1), . . .)) 

< A<j,A r u(t - 1)) + ^ /?(ui(t), - l))j . 

Using the notation 7 := A$At < 1 and continuing this estimate backward in time we get 

N 

p(u(t + 1), u(t)) < j 2 p(u(t - l),u(t - 2)) + 7 2 P(vi(l ~ 1), Vi(t - 2)) 

i=l 

2V 

+ 7 ^/>K(t),^(t-l)) 
1=1 

t AT 

< < 7 *p(«(l), «(0)) + Yl ^ E " fc + X )' " k )) 

k=l i=l 
t 

<7*p(«(l),«(0)) +2ATC , ^7 fc f~ fc+1 

fc=l 

< 7*p(«(l),u(0)) + 2NCt (max(7, 7 ))* +1 ^ 

at exponential rate (since max(7,7) < 1), which yields the result. □ 

Proof of Theorem 12. II Choose any element i G B. According to the definition of the map 7jg and to 

the property of the unidirectional interaction the value (7jgx)i depends only on X{ and the values at 
sites with smaller indices according to the ordering given by the function M. 

Our strategy is to apply Lemma fa . II consecutively (according to the enumeration TV) to all elements 
of the box B starting from the 'smallest' one (say i G B) having the number 1. In that case the 
corresponding value {T\ B x)i depends only on Xi and the values at sites outside of the box B given by 
boundary conditions. Thus Xi(t) is defined by the relation of type 1)3. ljl with A$ := A^. Hence by 
Lemma 13. II x; (t) converges to a certain limit time t — > 00. 
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Then we apply Lemma fa.ll to the site corresponding to the second smallest value of jV, etc. Observe 
that the exponential rate of convergence needed to apply the second assertion of Lemma 13.11 follows 
from the first one. 

It is important that in any given step of the procedure the values at sites with smaller indices 
(according to the enumeration Af) whose entries enter into the formula (|3.1jl and belong to the corre- 
sponding collection T do not depend on the value of the site under consideration. □ 

Corollary 3.2 Under the assumptions of Theorem \2.1\ the limit solution x depends continuously on 
the boundary conditions. 

Let us study how a limit solution depends on changes in boundary conditions. To simplify calcula- 
tions we consider here only the case of linear one-dimensional interactions given by (|2.4|) and identical 
one-dimensional maps Tj = T acting on a unit circle. We assume also that the boundary conditions 
are homogeneous, i.e. Xi = v Vi £ B c , and in some neighborhood U of the point uel the map T is 
linear and expanding, i.e. Tx = ax + b with \a\ > 1 for all x G U. Following the same argument as 
in the proof of Theorem 12. II we start with the site corresponding to the smallest positive value of the 
function J\f. Denote the value of the corresponding limit solution (existing by Theorem I2.1j) at this 
site by u. Then u satisfies the relation 

u = (1 — c)Tu + cTv . 
Substituting the local linear representation of T we get 

u = (1 — c)(au + b) + c(av + b) , 

or 

cav + b 

u = r . 

ca — (a — 1) 

Choose another homogeneous boundary condition v' € U. Then the corresponding limit solution can 
be written as 

, cav' + b 
u = . 

co — (a — 1) 

Therefore 

i /i ca , 

\U — U | = r • \V — V | 

ca — (a — 1) 

and the constant — 1 a ,\ > 1 since \a\ > 1 and (1 — c)a < A^At < 1. Note that we need to assume 
here that u, u' G U, so the neighborhood U should be large enough. 

Continuing this procedure for sites with larger TV-indices we see that the distance between limit 
solutions corresponding to the homogeneous boundary conditions v and v' grows exponentially with 
distance from the boundary of the box (at least until one of limit solutions leaves U). 

4 Proof of Theorem E2] 

First assume that the interaction in B is unidirectional but the graph G{B) has a cycle contrary to our 
claim. Then when we are moving along this directed cycle the enumeration function M{-) can only 
grow which contradicts to the existence of the cycle. 

Assume now that the directed graph G(B) has no cycles. We construct the enumeration function 
AA(-) as follows. We shall say that a vertex i € B is a starting one if there is no j £ B such that the 
edge (j, i) G G, and free if there is no j S B such that € G. 

Clearly due to the absence of cycles in G we deduce that the set of starting vertices is non empty. 

We start the enumeration procedure by choosing some starting vertex say i € B and setting 
M{i) := 1. 
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Assume now that we already have enumerated M vertices from B. Then we normalize the function 
M in such a way that for each k = 1, 2, . . . , M its fc-th largest value becomes equal to k (i.e. it is an 
exact numbering). 

Choose an arbitrary path i\,i2, ■ ■■ ,i n i n the graph G satisfying the properties that the function 
Af is not defined on all vertices of this path, the vertex %\ is starting or 3j G B : (j, ii) G G and Af(j) 
is well defined, and either (a) i n is free or (b) 3j' G B : (i n ,j') G G and Af(j') is well defined. 

If the condition (a) is satisfied we set 

N{i k ):=M + k, k = 1,2, . . . ,n, 
otherwise if the condition (b) is satisfied we set 

M(i k ):=M(j')-l + -^—, k = l,2,...,n. 
n + 1 

After this we normalize the function Af, choose a new path, etc., until no non-enumerated elements 
left in B. 

Our aim was to define the function Af on this path in such a way that the unidirectional interaction 
property should hold with respect to previously numbered vertices. It is straightforward to check that 
the enumeration obtained on each step of our construction does satisfy this property. □ 



5 Discussion of boundary conditions of different types 

The result of Theorem 12 . 1 1 mav suggest that at least in the case of homogeneous boundary conditions 
the mechanism leading to the LRA phenomenon is a simple propagation of the boundary conditions 
through the boundary into a box B. In general this is not the case and the limit solution x can be 
rather inhomogeneous in space. However this propagation phenomenon still may take place as we 
shall demonstrate in the following result. Assume that all local maps Tj are identical and denote this 
common local map by T. 

Lemma 5.1 Let the map T have a fixed point £ G X (i.e. T£ = £) and suppose that the boundary 
conditions are fixed at the value £, i.e. Xj(0) := £ Vi G B c . Then under the conditions of Theorem \2.1\ 
for any initial data x\jg(0) the solution x(t) will converge with time to the constant solution with 
{x)i = £ for all i G B. 

Proof. First, according to Theorem 12.11 there exists the only one limit solution x. On the other 
hand, by definition the interaction I preserves the "diagonal" and hence any sequence with identical 
coordinates is preserved by I. Therefore the vector with identical coordinates (x)i = £ for all i G B is 
invariant for both the map and the interaction and is, indeed, the only solution to our problem. □ 

It turns out that the solutions depend rather sensitively on the way how one defines the corre- 
sponding boundary conditions. 

Let us start with a very mild change in the definition of the interaction between the sites inside of 
the box B and the outer ones. Previously we used the convention that in order to calculate the value 
(T\jgx)i(t + 1) one first applies the map T at all sites and then applies the interaction 2 to the result 
of the previous operation. Contrary to this, assume that we apply the map only to the sites from the 
box B (but not to the outer sites) and denote the map describing the corresponding system by . 
From the point of view of Theorem 12 . 1 1 this change does not matter much and its claim remains valid 
for the modified map . Indeed, it is straightforward to check that the modified map with 
certain boundary conditions x'jgc is equivalent to the original map 7|g with the boundary conditions 
Xi G T~ 1 x' for all i G B c . Of course, we assume here that the preimages do exist, i.e. T~ 1 x' ^ 0. On 
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the other hand, the result about the propagation of homogeneous boundary conditions (Xemma l5.1|) is 
no longer valid. The problem is that the homogeneous vector x immediately becomes inhomogeneous 
under the action of . 

The boundary conditions considered so far can be called frozen because they do not change in 
time. We shall consider now two other types of boundary conditions. The first of them, called free, 
corresponds to the situation when after each application of the map the values at the sites outside of 
the box B do change under T, i.e. 

T^z) :=((( f ° f H if i G (51) 
' 1 \TiXi otherwise. 

Theorem 5.1 Let the interaction in a finite box B be unidirectional, the boundary conditions are free, 
for any i G B there be a path in G(B,I) from a site outside of the box B to i, and let A^At < 1. Then 
for any (initial) boundary conditions x# c (0) there exists a vector-function x(t) such that 

for any initial data x*b(0). 

Proof. Observe first that Lemma l3.1l has been proven for time dependent entries. Now Theorem 15. II 
follows essentially from the same argument as the one used in the proof of Theorem 12.11 This proves 
the existence of the vector- function x(t) which is the limit solution (depending on time now) to the 
problem with free boundary conditions. In fact the only difference to the frozen case is that the limit 
solution is no longer a constant. □ 



A completely different situation takes place if one considers periodic boundary conditions. Of 
course, in this case a box B must have a rectangular shape 2 (otherwise the periodic continuation is 
not well defined). Then the dynamics inside the box is completely determined by the initial data 
inside the box B. However, in distinction to the previous cases a limit solution might not be unique. 
Assume for simplicity that Tj = T for all i and the map T has two different periodic points £ and 
r) with non-intersecting trajectories. Consider two different initial data: a?j(0) := £, x'i(0) := rj for 
i £ B. Then the functions x{t) := T*£ and x'{t) : = T'ry describes two different time-periodic solutions 
to the problem having no common limit points. In fact, we even do not need any specific property of 
the interaction here and use only the invariance of the diagonal. 



6 Stochastic perturbations of local dynamics 

Conditions under which the existence of the LRA has been proven in Theorem 12.11 look very robust 
and certainly remain valid under small deterministic perturbations of local dynamics. 3 Our aim in 
this section is to show that the situation changes drastically if instead of deterministic perturbations 
one considers random ones. For a general introduction to random perturbations of chaotic dynamical 
systems we refer the reader to 0. However, since the interaction map I strictly speaking is not only 
non-contracting but non-hyperbolic as well, one cannot use general results here and everything has to 
be done "by hands". 

From now on we assume that all the assumptions of Theorem l2.1l are satisfied and hence the LRA 
phenomenon takes place for the unperturbed system. To simplify notation and to emphasize that 
the stochastic instability of the LRA is not caused by some artificial assumptions, throughout this 
section we assume that the local maps coincide (i.e. Tj = T) and the local phase space A is a unit 
(i-dimensional flat torus with the Euclidean metric on it. Note however that all results that we shall 
discuss remain valid when the maps Tj are different but satisfy to assumptions we shall make for T, 
and one can replace the unit torus by a smooth compact cf-dimensional manifold. 

2 or a more general polygon which one can fill in by reflections the entire plane. 
3 Namely, the vector x depends continuously on the local maps Tj. 
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Definition 6.1 A stochastic perturbation of a dynamical system (T,X) is a Markov chain with the 
same phase space X defined by the transition probabilities Q(x,A) to jump from a point x £ X to a 
Borel set A C X. By a stochastically perturbed system we mean a superposition of the original system 
and a perturbation, i.e. a new Markov chain defined by the transition probabilities Q(Tx, A). 

Definition 6.2 The amplitude of the stochastic perturbation is defined as the smallest number e such 
that Q(Tx, A) = whenever mm y£ A p(x, y) > £ for all x £ X. Therefore we shall use the notation Q £ 
to emphasize that the amplitude of the perturbation is s. We shall assume also that the transition 
probabilities are translationally invariant, i.e. for any x, z £ X and a Borel set A C X 

Q(x,A) = Q(x + z,A + z) . (6.1) 

Note that the phase space X is a torus here and hence the addition is well defined. 

To describe both deterministic and stochastic counterparts of the dynamics in the same terms it 
is useful to consider their action on the space of probabilistic measures M(X) on X. To distinguish 
between the pointwise action and the action on measures we shall use the upper index star for the 
latter (e.g. T* stands for the action on measures for the map T). These objects are defined as follows: 

T*fi(A) := up- 1 A), Qlix(A) := J Q £ (x, A)dfi(x) 

for any probabilistic measure [i £ A4(X) and any Borel set A C X. According to the previous definition 
an action of the stochastically perturbed system on measures is a product of these operators: Q*T*, 
while the action of the stochastically perturbed CML on measures can be written as 

f* :=l*Q*f* : M{X) -> M{X) , 

where Q* stands for the direct product of the operators Q* e corresponding to local stochastic pertur- 
bations. 

Definition 6.3 A solution jl £ £ M(X) to the equation T*jl = ft is called an invariant measure for 
the stochastically perturbed system. 

In distinction to our previous results which were concerned mainly with the stability type analysis 
we will need a number of assumptions to prove the instability against stochastic perturbations. Let us 
introduce the following notation. For \i € A4(X) denote by supp(/x) a support of the measure [i which 
defined as the intersection of all closed subsets of X of full /U-measure. We define also the collection 
of projection operators itk acting on the space M(X) where the index k £ Z rf corresponds to the k-th 
component of the direct product space X, i.e. 7Tfc/x stands for the projection of the measure jl £ M(X) 
to the k-th. component of the direct product of A4(X). 

We shall use the notation B r (x) for the ball of radius r centered at a point x £ B(X) and S x for 
the (5-measure at point x. To control a spread out of an initial measure fl £ M{X) under the action 
of the dynamics we introduce the following notion and a list of assumptions. 

Definition 6.4 By upper /lower spreads S±(fi) of a measure [i £ A4(X) we shall mean the radius of 
the smallest ball which contains support of the measure \x in the case of S + (//) , and the radius of the 
largest ball contained in the support of the measure /j, in the case of S-(/j,). 

The first of those functionals characterizes the maximal distance between points in the support of 
a measure, while the second one controls its concentration properties. 

We assume that there are constants At > Xt > 1 and Aq, a, b, a > such that 

Xrp{x,y) < p{Tx,Ty) < A T p(x,y) Vx,y : p(x,y) < a , (6.2) 

supp(Q £ *4) 5 B Xq£ (x) Vx £ X , (6.3) 
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S + (7r £ i*fI) > aS + (nifi) + b'Y]S + (ire>fI) (6.4) 

V 

for any i € B. The summation here and in the sequel is taken over all sites £' G such that 
(£' , €) € G. We assume also that the maximal degree of vertices of the graph G is bounded by K < oo. 

Theorem 6.1 Let the local maps, stochastic perturbations and the interactions satisfy the assumptions 
\6.ty) . \6.3\) . \6-4\ l and we consider a finite box B with certain frozen boundary conditions. Assume 
additionally that 

a\ T < 1, (a + b)\ T > 1. (6.5) 

Then the stochastically perturbed system in the box B admits an absolutely continuous invariant mea- 
sure fl e € M{X), and there exists a constant 7 > (which do neither depend on e > nor on B) such 
that 

S-(w e fle) >lL{t)e 
for each site I € B and small enough e > 0. 

In other words the lower spread of the invariant measure grows at least linearly with the 'distance' 
from the boundary L(£). Note the difference to a 'diffusive' growth which one might expect from 
random perturbations. 

Let now demonstrate that the assumptions (|6.2l - l63|) are natural. First note that the assumption 
Q6.2JI together with the assumption At > 1 imply only that the map T is locally expanding. Thus, 
e.g. the dyadic map Tx := 2x modi defined on a unit circle X considered as the interval [0, 1) with 
the metric p on it defined as the length of the shortest arc between points will satisfy them if we set 
Xt = At = 2 and a = 1/3. On the other hand, the assumption a Ay < 1 is necessary to be consistent 
with the LRA phenomenon in an unperturbed system (Theorem 12. 1(1 . 

To understand the reason behind the assumption (|6.4|) consider the simplest linear one-dimensional 
interaction defined by the relation (|2.4j) . for which we have a = 1 — c and b = c. Then Q6.4J) and the 
inequality (a + 6) At > 1 clearly hold true. 

The assumption Iffi.MJ) gives the information about the local stochastic perturbations of a point mass 
and clearly is satisfied in the case of independent local stochastic perturbations uniformly distributed 
in an e-neighborhood of a point with Aq = 1. 

Proof of Theorem 16. II Since the construction we shall apply is somewhat involved we start with the 
main scheme of the proof. Similar to the construction of the limit solution in the proof of Theorem l2.ll 
we shall construct the invariant measure \x e step by step starting from the first site according to the 
enumeration J\f and adding each time only one site according to this enumeration. 4 Using that the 
dynamics on already considered sites do not depend on the dynamics of the new one (but not vice versa) 
we may (and will) consider the interaction with those sites as an additional stochastic perturbation. 
On this stage our main goal is to show that the action on measures of the stochastically perturbed 
system is a contraction in a suitable metric in the space of measures. To this end we shall show that 
the expansion rate of the operator T* is proportional to At, while for stochastic perturbation Q* it 
is close to 1. The key point is to show that the expansion rate of projections of the operator X* are 
proportional to A^, which together with the assumption A^A^ < 1 guarantees the contraction of the 

operator T* and hence the convergence in time to the unique invariant measure. On the next step of 
the construction we shall obtain an upper bound for the upper spread of the projection of the invariant 
measure fl e to the £-th subsystem as a function of the 'distance' to the boundary L(£) and then using 
it we shall get the lower bound for the lower spread of this projection. 

We introduce the so called Hutchinson metrics in the space of probabilistic measures M(X) defined 
by the relation 

p*(fi,u) := sup ( / (pdfi — / ipdv) , 
K v <\ J J 

4 In fact, the limiting solution constructed in the proof of Theorem 12.11 describes the support of the invariant measure 
/2o of the unperturbed system. 
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where the supremum is taken over all continuous functions 99 : X — > R 1 satisfying the assumption that 
their Lipschitz constants A^, (defined by the relation do not exceed 1. 

Lemma 6.1 Under the assumptions of Theorem \6.1\ for any pair of measures p,r] € A4(X) we have 

p*(T*p,T*v)<A%p*(p,v) , (6.6) 

p*{Q* e p,Q»<P*M , (6.7) 
p*(Tr e i*p,TT e i*v) < A lP *(p,u) , (6.8) 

where ft, V G M.{X) are arbitrary measures satisfying the assumption that -Kip = p and -n^V = v. 
Proof. By definition 

p*(T*p,T*v) = sup(/ ipdT*p- [ <pdT*v) 
k v <\ J J 

= sup (/(/?o Tdp — ip o Tdv) . 
A<p<l J J 

Since p(Tx, Ty) < Atp(x, y) for any pair of points x, y S X we have 

\tp o T(x) -<poT(y)\ < A v p(Tx,Ty) < A (p A T p(x,y) . 

Thus setting ip := A^ip o T we get Av, < and hence 

p*(T*p,T*u) < A T sup ( / ipdp - / tpdu) = A T p*(p,u) , 

A^<1 J J 

which proves inequality (|6.6|) . 

Let us deal now with stochastic perturbations. 



Aip 

Using the translation invariance of Q E (-, •) we get 



p*(Q* e p,Q*v) = sup ( / ipdQ*p- / ipdQ* e v) 
A<fi<l J J 

sup ( / Q e <pdp - \ Q £ ipdu) . 

V(Z><1 J J 



\Q £ (p(x) - Q £ ip(y)\ =| / (p(z)Q £ (x,dz) - J (p(z)Q £ (y,dz)\ 

ip(z)Q £ (x, dz) - J (p(z - x + y)Q £ (y, dz-x + y)\ 
(p(z)Q E (x,dz) - J (p(z - x + y)Q £ (x,dz)\ 
< I \tp(z) -(p(z-x + y)\ -Q £ (x,dz) < A^p{x,y) . 



Thus setting if) := Q £ (p we obtain 

p*(Q* e p,Q* £ v) < sup ( / ipdp- / ipdv) = p*{p,v) , 
Ai/><1 J J 

which proves the inequality (|6.7() . 

To prove the remaining inequality (|6.8() we introduce the following notation. For a given x € X 
let xi Xl =x) £ X mean a vector whose l-th coordinate coincides with x, and for p£ := imp £ A4(X) 
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denote by pjr the projection of the measure p E M(X) to all subsystems except the £-th one (i.e. we 
integrate over the ^-th coordinate). Then 

(p(x)dTr e i*fx(x) = J (J ip((lx {xe=x) )) dfii(x {xe=x) )^ dpg(x). 

Consider a function 

ip(x) := AZ 1 / ip((lx {xi=x) )) dpf(x {xi=x) ). 



Then 



A f ■ \^{x 



) - 1>(y)\ < J \<p((ix(x e =x))) - <p((ix(x e = y )))\ dvi(x (xe=x) ), 



since pf{x^ Xl=x -^) = pf{x^ Xe=y ^) because the measure pf does not depend on xg. 
On the other hand, by the inequality 12.21 we have 

p(ix (xe=x) , lx {xi=y) ) < kfp{x,y) , 

which implies that 

Therefore using the same argument as in the previous estimates we obtain the inequality (|6.8|) . □ 



For a given < £ B the set T(£) consists of a finite number of elements connected to £ (directly 
or indirectly 5 ). Consider a subsystem consisted only of the sites belonging to T(£). The restricted 
dynamics of this subsystem does not depend on the dynamics of the other parts of the system. 

Collecting inequalities (|6.6|) . (|6.7|) . (|6.8|) we see that the operator corresponding to the stochastically 
perturbed system is strictly contracting in the p* metric, which (due to the compactness of this space) 
immediately implies that for any initial measure p G A4(X) we have: 

Going through all sites of the box B we construct the invariant measure p e for the entire stochas- 
tically perturbed system restricted to the box B with the given frozen boundary conditions. 

Now we are ready to obtain estimates from above for the upper spread S + (irgp e ). It is straightfor- 
ward to check that under (|6.2j) for any measure p S A4(X) whose upper spread satisfies the inequality 
S + (p) < a /At we have 

S + (T*p) < A T S + {p) . (6.9) 
The relation ()2.2|) implies that if S+(irgp), Sj r {'K^p) < a j At then 

I' 

Since the amplitude of the stochastic perturbations is bounded from above by e > 0, one has 

S + (Q* £ p)<S + (p) + e . (6.11) 

The dynamics at the site £' does not depend on the dynamics at the site £ (but not vice versa) therefore 
we may choose an initial measure p E Ai(X) satisfying the following conditions: 

S+fajZ) < a/(KA T ), TT k p = n k p £ Vfc : M(k) < M{£). 
5 through some other elements of B 
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Then by the inequalities (J6.91 16.101 I6.11|) we have: 

S + (n e f £ *fi) < A f (A T S + (^/2) + e) + A f ^(A T S+(^/2) + e) . (6.12) 

f 

Denote q := A^At, p := Applying the inequality l)6.12j) t times and using that S+(tt£>T* jl) = 

S+faijl) we get 

S+faftYfi < qtS+fafl + (q^S+fafi +pe) ■ 

Now passing to the limit as t — > oo and using that g = A^At < 1 (by the condition of Theorem 12. 1[) 
we obtain the desired estimate from above: 

q \ - , _ . 



Thus the upper spread of the £-th projection of the invariant measure cannot exceed a constant in 
power L(l) multiplied by e. Observe that for each t > 1 the inequality Z)(7r^(7^*)*/2) < cr holds 
automatically since g < 1. 

Now whence we have established the control over the upper spread of ji e we will obtain an estimate 
from below for the lower spread of this measure. 

Under (|6.2() for any measure [i € M(X) with S-(fi) < a/ At we have: 

S_(T» > X T S-(fi) . (6.13) 

The relation ()6.3|) implies 

S-(Q* £f i)>S^) + \ Q e. (6.14) 

Making use of the inequalities ()6.13j) . (|6,14j) and ()6.4|) for any measure \i G A^(X) with S-{\i) < a/ At 
and a measure /I € .M(X) with irijl = \i we get: 

S-{^f*p) > a(\ T S-(Trtfl) + A Q e) + 6(A r + Age) . (6.15) 

Denote a := aAr < 1, := 6Ay > 1 — a (by the assumption (|fi.5j) ) and k := (a + 6)Ag. Now, 
applying (|6.15|) to (7^.*)* we have 

S-inii^Yp) > afS-(nP) + (/3 J>-(w) + K£ ) 4" n 



Passing to the limit and using that /?/(! — a) > 1 we obtain the final estimate 



a 



Theorem 16. II follows by setting 7 := ke/(1 — a). □ 



7 Conclusion 

We have shown that the absence of loops in directed connectivity graphs of networks of dynamical ele- 
ments results in appearance of a long range action provided that interactions between local dynamical 
systems (elements of a network) are sufficiently strong. If local systems have a regular (contractive) 
dynamics than the LRA appears already for arbitrarily weak interactions (conditions of Theorem 12.11 
clearly hold in this case). 
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One expects that this restriction to the network topology is by no means a necessary one. In 
principle, we see two ways to overcome this problem. The first of them is to use some rather general 
ergodic-theoretical results (see e.g. 0) to recode interaction potentials with one-sided dependence to 
the ones with the two-sided dependence. It is known that if those potentials are homological to each 
other than they generate the same Gibbs measures (corresponding to our invariant measures). The 
problem here is that it is hard to describe the set of constructed potentials corresponding to general 
non- unidirectional interactions. 

Another possibility is to consider non-unidirectional interactions as small perturbations to the 
unidirectional ones. This approach seems to be more promising but one inevitably needs to make a 
number of additional assumptions about the unperturbed system to ensure that the unidirectional 
interactions are stronger than the perturbations. 

There is yet another possibility related to the assumption that the interaction in the network is 
close to a mean-field model, i.e. a large part of the network is interacting with each its element. The 
main disadvantage of this assumption is that it is not completely clear how to extend this property 
for infinite networks. 

It is worth mentioning that instead of the lattice Z rf one can consider a countable ordered collection 
of vertices. The short range interactions in this context would mean that subsystems corresponding 
to vertices interact only if the distance between them according to the above mentioned order does 
not exceed a certain threshold. Clearly all our results remain valid in this setting. 

The problems similar to the ones discussed in this paper were studied in the context of the synchro- 
nization theory. More specifically, we are talking here about a complete synchronization of dynamical 
systems jS]. The complete synchronization means that asymptotically, as time goes to infinity, the 
differences in the behavior of all elements in the network disappear. In particular, in |1U| it has been 
shown that in the case of the network described by the ordinary differential equations whose directed 
connectivity graph is a tree the synchronization phenomenon takes place under strong enough inter- 
actions, while under the assumption that the connectivity graph is "strongly connected" (i.e. the 
interaction is close to a mean-field model) this has been proven in Tj (see also references in these 
papers and the recent monograph [H]). In these publications, in distinction to our approach, a syn- 
chronization phenomenon is studied in closed finite networks consisting of identical elements and the 
synchronized system may demonstrate chaotic behavior. 
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